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1. Introduction 

Understanding the deconfinement phase transition is one of the major challenges of particle 
physics. In quenched QCD reliable results are obtained within the lattice approach. This approach 
fails, however, at large baryon density due to the notorious fermion sign problem. Therefore alter- 
native non-perturbative approaches to continuum QCD are desirable. In recent years a variational 
approach to Yang-Mills theory in Coulomb gauge was developed [|l|], which has provided a decent 
description of the infrared sector of the theory [|[ |], ^ |[ |], ^. Recently this approach was ex- 
tended to finite temperature [^] and also to full QCD [^. In this talk I will report on the calculation 
of the effective potential of the confinement order parameter within the Hamiltonian approach [p^]. 

In quantum field theory the temperature T is most easily introduced by compactifying the Eu- 
clidean time and interpreting the length L of the compactified time interval as inverse temperature. 
In finite temperature Yang-Mills theory the order parameter of confinement is the expectation value 
of the Polyakov loop 

P[Ao] = itr^e--^«''^^"^''(^""^) . (1.1) 

The quantity {P[Ao](x)) ~ exp [—Foo{x)L] is related to the free energy of a (infinitely heavy) quark 
at spatial position x. In the confined phase this quantity vanishes by center symmetry while it is 
non-zero in the deconfined phase, where center symmetry is broken. In continuum Yang-Mills 
theory the Polyakov loop is most easily calculated in Polyakov gauge ^o^o = , Aq = diagonal. 
In the fundamental modular region < AqL/I < n the Polyakov loop P[Ao\ is a unique function 
of the field Aq, which, for SU(2), is given by P[Ao] = cos(AoL/2). As a consequence of this 
relation and of Jenssen's inequality one can use instead of (/'[Aq]) alternatively P[(Ao)] or (Aq) as 



order parameter of confinement, Refs. [|nj, |12|]. The order parameter of confinement can be most 
easily obtained by calculating the effective potential e[aQ] of a temporal background field oq chosen 
in the Polyakov gauge and by calculating the Polyakov line ( |1.1[ ) from the field configuration ao 
which minimizes e[ao], i.e. (P[Ao]) ~ ^[^o]- The effective potential e[ao] was first calculated in 
Refs. [13, 14] in 1-loop perturbation theory and is shown in Fig. [T[ This potential is minimal for 



a vanishing field and the order parameter accordingly yields P[dQ = 0] = 1 , which indicates the 
deconfining phase. The aim of the present work is to give a non-perturbative evaluation of e[ao] 
[10] in the Hamilton approach to Yang-Mills theory 

It is obvious that the effective potential of (Aq) cannot be straightforwardly evaluated in the 
Hamiltonian approach since the letter assumes Weyl gauge Aq = 0. However, we can exploit 0(4) 
invariance of Euclidean quantum field theory and compactify instead of the time one spatial axis 
(for example the XT,-a.xis) to a circle and interpret the length L of the compactified dimension as 
inverse temperature. Therefore we will consider in the following Yang-Mills theory at a finite 
compactified length L in a constant color diagonal background field as and calculate the effective 
potential ^[03]. In the Hamiltonian approach the effective potential e[a] of a spatial background field 
a is given by the minimum of the energy density {H) /V calculated under the constraint (A) = a. 
This minimal property of the effective potential calls for a variational calculation. 

2. Hamilton approach in background gauge 

In the presence of an external constant background field a the Hamiltonian approach can be 
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most conveniently fomiulated in the background gauge 



d,A 



0, d = d +a 



(2.1) 



where all fields are taken in the adjoint representation. This gauge allows for an explicit resolution 
of Gauss' law, which results in the gauge fixed Hamiltonian 

1 



H = -jd\ [J^'U{x)Ja ■ h{x) +B\x)j +Hc , 

where ft = —i5 / 5A is the momentum operator of the gauge fixed field and 

JA=Det(-D-d) , D = d+A 



is the Faddeev-Popov determinant. Furthermore, 

,2 



//c = |-/ d\d'yJ^'p"{x)JAF"'{x,y)p'{y) 



(2.2) 



(2.3) 



(2.4) 



is the analogue of the so-called Coulomb term which results from the kinetic term of the "longitu- 
dinal" part of the momentum operator. Here 



A-a] n 



-DU 



is the color charge density of the gluons, which interacts through the kernel 



-Dd 



-d-d 



-Dd 



(2.5) 



(2.6) 



For a vanishing background field a = the gauge ( [2. Ij ) reduces to the ordinary Coulomb gauge and 
the Hamiltonian H QJl) becomes the familiar Yang-Mills Hamiltonian in Coulomb gauge, [p^]. 
Furthermore, as was shown in Ref. [||] the Coulomb term Hq (2.4) is negligible in the gluon sector. 
Therefore we will ignore this term in the following. 

We are interested here in the energy density in the state [^] minimizing 



{H)^:={yifa\HWa) 



(2.7) 



under the constraint {A) a = a. For this purpose we perform a variational calculation with the trial 
wave functional 

¥a[A] = Ja''^W[A - a] , \ff[A] = ^e-5/A«A ^ (2.8) 

which already fulfills the constraint {A)a = a. For a = this ansatz reduces to the trial wave 
functional used in Coulomb gauge Proceeding as in the variational approach in Coulomb gauge 
[|l|], from {H)a —5- min one derives a set of coupled equations for the gluon and ghost propagators 



, . 1 
{AA)a=0 =-(0 



-1 



G 



(D + a)d 



(2.9) 



a=Q 



Using the same approximation as in Ref. ^ in Coulomb gauge, i.e. restricting to two loops in 
the energy, while neglecting He and also the tadpole arising from the non-Abelian part of the 
magnetic energy, one finds from the minimization of {H)a the gap equation 



CO 



-d-d + x^, 



(2.10) 
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where' 



is the ghost loop (referred to as "curvature") with Fq and F being the bare and full ghost-gluon 
vertex. The gap equation ( 2.10D has to be solved together with the Dyson-Schwinger equation 



(DSE) for the ghost propagator 

= -d- d- Fo(l)GF(2)^(2, 1) . (2.12) 
Due to the presence of the background field these equations have a non-trivial color structure. 



Fortunately, due to the choice of the background gauge ([2. 1[), the background field enters these 
equations only in form of the covariant derivative d = d + a. 

We are interested in the effective potential of a background field in the Cartan algebra, which 
for SU(2) has the form aT^. This field becomes diagonal in the Cartan basis defined by the eigen- 
vectors of the generators of the Cartan subgroup Tsja) = — ia|a) , a = 0,ibl. By analyzing our 
equations of motion, i.e. the gap equation and the ghost Dyson-Schwinger equation, one finds that 
these equations have solutions, which are diagonal in the Cartan basis 

^'^^p) = , G^"(75) = 5^"G^(^) . (2.13) 

This is not surprising since the source of the non-trivial color structure is the background field and 
if this field is chosen to be diagonal in color space the same should be true for all propagators. In 
addition, one can show that the propagators in the presence of the background field in background 
gauge are related to the propagators in Coulomb gauge (i.e. in the absence of the background field), 

^{p), G{p), by 

^-(p) = ^(p-), G-(p) = G(p-), (2.14) 

where 

p''=p-oa (2.15) 

is the momentum shifted by the background field. In this way the results of the variational calcu- 
lation in Coulomb gauge (in the absence of the background field) are sufficient to determine the 
propagators in the presence of the background field in background gauge. 

Lattice calculation [ p^ ] of the gluon propagator in Coulomb gauge show that the gluon energy 



can be nicely fitted by Gribov's formula [17] 



(0{p) = \ p^ + M^p^. (2.16) 



A full self-consistent solution of the gap equation ( |2.10| ) and the ghost DSE ( 2.12| ) reveals that 



(o{p) contains in addition sub-leading UV-logs, which on the lattice are found to be small. Using 
Gribov's formula ( |2.16[ ) for a){p) and solving the gap equation ( |2.10| ) for x{p) yields 



Xip)=M^/\p\, (2.17) 

which is indeed the correct IR-behavior obtained in a full solution [||] of the coupled ghost DSE 
and gap equation but which misses the sub-leading UV-logs. 

'We use here the compact notation = A"^' (xi). For Lorentz scalars hke the ghost, the index "1" stands for the 
color index ai and the spatial position .vj . Repeated indices are summed/integrated over. 
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Figure 1: The Weiss potential eyv multiplied by 
L"* as a function of the dimensionless field x — 
aL/{2%). 



Figure 2: The infrared potential eiR multiplied 
by L!^ for different L ' . 



3. The effective potential 



Compactifying the 3-axis to a circle with circumference L and choosing the background field 
along the compactified dimension a = ae^, the shifted momentum ( 2.15D becomes 



p'^ =p± + iPn-<ya)e3, Pn = 27ln/L, 



(3.1) 



where p± is the projection of p into the 1-2-plane and p„ is the Matsubara frequency. In the 
Hamiltonian approach the effective potential of the constant background field is given by the energy 
density in the state minimizing {H)^ under the constraint (A)^, = a |18]. Using the gap equation 
2. 10| one finds for the energy density per transversal degree of freedom in the present approximation 



[10] 



(2n) 

By shifting the summation index n one verifies the periodicity 

e{a + 2K/L,L) = e{a,L) , 



(3.2) 



(3.3) 



which is a necessary property for the effective potential of the confinement order parameter by 



center symmetry. Neglecting xip) Eq. (3.2) gives the energy of a non-interacting Bose gas with 
single-particle energy (o{p). This quasi -particle picture is a consequence of the Gaussian ansatz 



(2.8) for the wave functional. The quasi-particle energy (o{p) is, however, highly non-perturbative. 



see for example Eq. (|2.16| ). The curvature x{p) Eq. ( p.2|) arises from the Faddeev-Popov deter- 



minant in the kinetic part of the Yang-Mills Hamiltonian (2.2) 



In certain limiting cases and for < fl'L/27r < 1 the energy density ( p^ can be calculated 
analytically. Neglecting x{p) assuming the perturbative expression for the gluon energy 
(o{p) = \p\ one finds from (3^) the Weiss potential originally obtained in [ p^ ] 



AK^faLVfaL 
^uv(a,L) = -^(^-J (^--1 



(3.4) 



5 



Hamiltonian Approach to QCD: The effective potential of the Polyakov loop 



Hugo Reinhardt 




shown in Fig. |l|. Neglecting x{p) ^rid using the infrared expression for the gluon energy (o{p) 
M^/\p\ (see Eq. ( |2.16| )), one obtains |[T^] 



eiR{a,L) = 2jj 



aL\- 
Ik) 



aL 
In 



(3.5) 



shown in Fig. g. This expression drastically differs from the Weiss potential (3.4): While e\jY{a,L) 
is minimal for a = 0, the minimum of eiR(a,L) occurs at a = n/L corresponding to a center sym- 
metric ground state. Accordingly ei]Y{a,L) yields for the Polyakov loop (P) = P[Ao = 0] = 1 while 
eiR{a,L) yields (P) = P[Ao = %/L]= 0. 

The deconfinement phase transition results from the interplay between the confining infrared 



potential, Eq. ( |3.5[ ) and the deconfining UV-potential, Eq. (|3.4|). To illustrate this let us approximate 
the gluon energy co{p) ( |2.16[ ) by 

co{p)^\p\+M'-/\p\. (3.6) 



This expression agrees with the Gribov formula ( [2. 16| ) in both, the IR and UV but deviates from 
it in the mid-momentum regime, which influences the deconfinement phase transition. With (o{p) 



given by Eq. (g^ and with % (p) 

e{a,L) = eiR{a,L) + e\jY{a,L 



the energy density ( [3.2| ) becomes 
4n^ 



-J 



faL\ 



f{x) =x^{x- \f + cx{x- 1), 



3L4-' ^271 J ' 
_ 3M2l2 



(3.7) 



For small temperatures e\R{a,L) dominates and the system is in the confined phase. As L^^ 
increases the center symmetric minimum at x = 1/2 eventually turns into a maximum and the 
system undergoes the deconfinement phase transition, see Fig. ^ In the deconfined phase f{x) has 
two degenerate minima and, starting in the deconfined phase, the phase transition occurs when the 
three roots of f'{x) degenerate. This occurs for c = 1 /2, i.e. for a critical temperature 



\/3M/7r. 



(3.8) 
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Figure 5: The numerically evaluated effective 
potential. The critical temperature is ~ 269 
MeV. 
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Figure 6: The Polyakov loop {P[a]) evaluated at 
the minimum a = amin of the full effective poten- 
tial (figure (a)) as a function of T/Tc- 



With the lattice result M = 880 MeV this corresponds to a critical temperature of ~ 485 MeV, 
which is much too high. This value is only slightly reduced to = 432 MeV when the correct 
Gribov formula ( 2.16 ) is used instead of the approximation (3^). The reason why the transition 
temperature comes out too high is the neglect of the ghost loop x{p) = 0. This can be seen from 
Fig. ^ where the gluon energy 0){p) and the curvature x{p) obtained from the variational calcula- 
tion in Coulomb gauge are shown. In the deep infrared the gluon energy 0){p) and the curvature 
X{p) agree, while in the ultraviolet x{p) has the opposite sign of 0){p) and is only suppressed by 
a log compared to 0){p). Therefore, neglecting the ghost loop enhances the infrared contribution 
to the potential, which favors confinement, while at the same time it reduces the ultraviolet con- 
tribution, which favors deconfinement. Both effects add coherently and push the deconfinement 
transition to a higher temperature. 

In a full numerical evaluation of the effective potential (3^) using for 0){p) and x{p) the 
numerical solution of the variational approach in Coulomb gauge obtained in Ref. [||], one finds 
the effective potential shown in Fig. From this potential one extracts a critical temperature for 
the deconfinement phase transition of Tc ~ 269 MeV, which is close to the lattice predictions of 
Tc = 290 MeV. This value is also close to the range of critical temperatures Tc = 275 . . . 290 MeV 
obtained in Ref. [||] from the grand canonical ensemble of Yang-Mills theory in Coulomb gauge. 
Let us also mention that if one uses for co{p) the Gribov formula ( 2.16D and in accord with the gap 
equation (2.1C) for x{p) its infrared expression (2.17) one finds a critical temperature of Tc ~ 267 
MeV, which is only slightly smaller than the value obtained with the full numerical solution for 
(o{p) and x{p)- This shows that it is the infrared part of the curvature (neglected in Eq. (3.7)), 
which is crucial for the critical temperature. In view of the ghost dominance in the IR this is not 
surprising. Fig. ^ shows the Polyakov loop P[a] calculated from the minimum amin of the potential 
(3^). At the phase-transition /'[amin] rapidly changes from P = to P = 1, which is typical for a 
second order phase transition. 

In the present approach the deconfinement phase transition is entirely determined by the zero- 
temperature propagators, which are defined as vacuum expectation values. Consequently, the finite- 
temperature behavior of the theory and, in particular, the dynamics of the deconfinement phase 
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transition must be fully encoded in the vacuum wave functional. The results obtained above are 
encouraging for an extension of the present approach to full QCD at finite temperature and baryon 
density. 
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